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We study, in the ensemble of constant particle number, processes in which a cluster of particles 
is annihilated and particles are created catalytically in active sites. In this ensemble, particles 
belonging to a cluster of £ particles jump to £ distinct active sites. As examples of our prescription, 
we analyze numerically three nonequilibrium systems that annihilate cluster of particles that are 
identified as conserved versions of the pair annihilation contact model, triplet annihilation contact 
model and pair contact process. We show also how to set up the constant particle number ensemble 
from the constant rate ensemble. 
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I. INTRODUCTION 

The use of distinct ensembles in equilibrium statisti- 
cal mechanics, in which systems are described by a given 
Hamiltonian, is a well established concept There is 
a standard procedure for passing from a given ensemble 
to another. For nonequilibrium systems on the other 
hand, this procedure can not be used, since they are not 

] described by a Hamiltonian and therefore, their proba- 
bility distribution are not known a priori. In most cases 
nonequilibrium systems are defined in a constant rate 
ensemble, (which we call ordinary version), but the pos- 
sibility of using another ensemble was put forward by Ziff 
and Brosillow 3] when they used a constant coverage en- 
semble in their study of an irreversible surface-reaction 
model. Subsequently, Tome and de Oliveira [3| intro- 
duced the contact process in the ensemble of constant 

' number of particles. 

In the conserved contact process (CCP) [1], a particle 

I chosen at random leaves its place and jumps to one of the 
many active sites of the lattice. The CCP displays prop- 
erties that in the thermodynamic limit are identical to 
the ordinary contact process. Hilhorst and van Wijland 
[H provided a proof of the equivalence between the con- 
stant rate and the constant particle number ensembles 
for the contact process. Later, de Oliveira extended 
the proof for any reaction process that annihilate one 
particle 0, 0, 0- Sometimes the use of the conserved 

I version is more appropriate than the ordinary version as 

, for example in the study of a first-order transition. This 
advantage has been exploited by Ziff and Brosillow 3, 
Loscar and Albano ^ that studied hysteretic effects in 
a model that describes the CO+NO reaction, and more 
recently by Fiore and de Oliveira {8]. 

Here, we analyze in the ensemble of constant particle 
number one-dimensional models whose ordinary versions 
have been previously studied. These models are: the 
pair-annihilation contact model (PAM) 11 Oil , the triplet- 
annihilation contact model (TAM) (Tol . Illl |. and the pair 
contact process (PCP) [T2|. ITsl. [T4| . In these models, a 
cluster of particles is spontaneously annihilated and par- 



ticles are catalytically created in active sites. Active sites 
are empty sites surrounded by a neighborhood of parti- 
cles. 



II. TRANSITION RATES 

Consider a site i of regular lattice. To each site i of the 
lattice we attach an occupation variable tji which takes 
the values or 1 according whether the site i is empty or 
occupied by a particle. In the constant rate ensemble, the 
usual process is composed of creation of a single particle 
(0 —>■ 1) with transition rate = kcwf, and annihilation 
of a cluster of £ particles in a row (111...1— »000...0) with 
transition rate = kaUjf. The transition rate wf is 
the probability per unit time of creating a particle at the 
site i. The transition rate wf is the probability per unit 
time of annihilating a particle at the site i. The total 
transition rate for these two reaction processes is given 
by 

Wi = Wi + Wi = kaUJi + kcLof. (1) 

The quantities ujf and will be defined according to the 
specific model and they depending on the local configu- 
ration of particles. The quantities ka and kc, which we 
call strengths of the annihilation and creation process, 
are parameters that gives the weight of each subprocess. 
Diffusion of particles consists of a particle hopping to its 
nearest neighbor. In diffusive models, the diffusion of 
particles and the reaction process occur with probability 
D' and 1 — D' , respectively. 

In the ensemble of constant particle number, the jump- 
ing process and the diffusion of particles are chosen with 
probabilities 1 — D and D, respectively. In the hop- 
ping step, an occupied site rji = 1 and its nearest next 
neighbor r]j = interchange their occupation variables, 
whereas in the jumping step, £ adjacent particles leave 
their places and arrive at £ active sites. Thus, in the 
ensemble of constant particle number, both the creation 
and annihilation of particles are replaced by £ jumping 
process. 
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III. EQUIVALENCE OF ENSEMBLES 

In the construction of conserved models we have to be 
concerned only with the creation and annihilation pro- 
cesses since the diffusion process already conserves the 
particle number. In the case of processes that annihilate 
a pair of particles (PAM and PCP), the jumping process 
occurs with rate wfcj^w^. In the case of annihilation of 
a triplet of neighboring particles (TAM), it occurs with 
rate w>j^w^cj^. 

To demonstrate that this conserved dynamics is equiv- 
alent to the ordinary dynamics described by Ekj. ([!]) we 
follow the same reasoning given by de Oliveira ^6] . Let us 
consider the specific case of models that annihilate two 
particles. In this case, two particles, belonging to the 
neighborhood of site i, chosen at random, jump to two 
distinct active sites, j and k, also chosen at random. The 
jumping process occurs with rate = ujftOjLLj'^/L'^ 

where L is the number of sites of the lattice. Let us eval- 



uate the total rate J2 



j,k 



in which particles leave 



the neighborhood of site i. In the thermodynamic limit, 
the sums J2j ^j/^ ^^-'^ J2k^k/-^ approach, by the law of 
large numbers, the averages (ujj) and (ij-'j:), respectively. 
Therefore we have J^jk^iJ^k = {'-^j){'-^k)^i- Similarly, 
the total rate at which a particle arrives at the site j is 
{u!f){Lu^)u!j and the total rate at which a particle gets to 
the site k is {lo°;){u)'^)uj1. 

Comparing these results with the rate ([I]), we see 
that kc and ka should be proportional to 2{ujf){ujl) and 
(a;|)(a;£), respectively. Defining a as the ratio ka/kc, we 
can write the following relation 



(2) 



valid for any process that annihilates a pair of particles 
in the ensemble of constant particle number. 

When we add a diffusive step, the ensembles are also 
equivalent, but the diffusion rate D' used in the constant 
rate ensemble will not have, in general, the same value 
of diffusion rate D that is used in the constant particle 
number ensemble. A relation between the rates D' and 
D is given by 



D' 



D 



P P 



(3) 



where the factor 1/ comes from the ratio between the 
occurrences of diffusive process and the jumping process. 

Generalizing for an arbitrary cluster of £ particles, the 
Eqs. ^ and Q become 



Ml 



and 



l-D' 1- D 



D' 



D 



(4) 



(5) 



respectively. In the particular case of £ = 1, we have 
D = D because (wf) = {-qi) = p, as it has already been 
obtained and used in the simulations of diffusive contact 
processes @, 



IV. PAIR-ANNIHILATION CONTACT MODEL 

A. Constant rate ensemble 

In the ordinary PAM, the creation of particles is cat- 
alytic and a pair of particles is annihilated spontaneously. 
It is represented by the chemical reactions: 



+ A ^ ^ + A, 



(6) 



(7) 



describing the catalytic creation and annihilation of a 
pair of particles, respectively. The rates cjf and cjf in 
Eq. ([!]) are given by 



(8) 



where the summation is over the z nearest neighbor sites 
and 



(9) 



The strengths of the creation and annihilation processes 
are given by fcc = 1 and ka — a. 

For values of a > Ofc, the system is constrained into the 
absorbing state (without particles), whereas for a < ac 
we have an active state in which particles are created and 
pairs of neighboring particles are annihilated. A continu- 
ous phase transition between these two regimes occurs at 
a = ac = 0.18622(3) pX?|. Close to the critical point, the 
order parameter (here the density of particles p) follows 
a power law 



p - (ac - a)^. 



(10) 



where (i = 0.2765(1). The PAM is found to belong to 
the directed percolation universality class (DP) 0. 



B. Constant particle number ensemble 

The Monte Carlo simulation of the conserved PAM is 
performed as follows. The jumping process and the dif- 
fusive process are chosen with probabilities 1 — D and D, 
respectively. If the jumping process is chosen, then a pair 
of neighboring particles is chosen at random. Next, we 
choose two empty sites each one surrounded by at least 
one particle. The two particles belonging to the selected 
pair jump to these two active sites. If the diffusive pro- 
cess is chosen, then a particle chosen randomly is moved 
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FIG. 1: Values of a versus the number of particles n for the 
conserved PAM in the absence of diffusion for an infinite sys- 
tem. The line corresponds to the linear extrapolation using 
Eq. (HlJ. 

to one of its neighboring site, provided it is empty. In 
this ensemble, the rate a is evaluated by using formula 
([2]), where the quantities wf and are given by Eqs. ([8]) 
and ([5]). To determine the critical rate ac we have sim- 
ulated the conserved PAM in the subcritical regime. In 
this regime we have an infinite lattice and a finite particle 
number n. The critical value ac is obtained by assuming 
the asymptotic relation 

1 / ^ 

a — ac ~ — ■ (11) 

A linear extrapolation of a versus gives Uc, since 
when 71 ^ CXI we have a — > ac- The behavior of a versus 
for the conserved PAM in the subcritical regime is 
shown in Fig. [TJ In the limit of n oo, we obtain 
ac = 0.18624(7), in excellent agreement with the value 
ac = 0.18622(3) obtained from its ordinary version p^ . 

According to Vicsek [l^, in the critical point, we have 
a formation of fractal clusters. To calculate the fractal 
dimension for a fixed D, we have measured the maxi- 
mum distance R between two particles of the system as 
a function of the particle number n. 

Following Broker and Grassberger [I^ , we assume the 
following asymptotic behavior 

Rr^n^/'^", (12) 

where dp is the fractal dimension. The fractal dimen- 
sion is related to the survival probability exponent S, the 
mean number of particles exponent 77, and the dynamic 
exponent z by dp = 2(S + rj)/z 2]. In Fig. [2] we show 
a log-log plot of R versus n for the conserved PAM for 
some values of the diffusion rate D. The values of l/dp 
are consistent with the value 1.338(6), obtained for the 

CCP [zi]. 



FIG. 2: Log-log plot of the size system R, in the subcritical 
regime, as function of n where n is the particle number n for 
the conserved PAM and the conserved TAM, and the number 
of pairs of neighboring particles rip for the conserved PCP. 
For comparison, we show a straight line with slope 1.338. 
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FIG. 3: Phase diagram for the conserved PAM in D versus a 
space. The supercritical and subcritical regimes are separated 
by a critical line. 

Using Eq. (fTTj) we have built the phase diagram shown 
in Fig. [3]for several values of diffusion rate D. Increasing 
the diffusion rate, we expect an increase in the value of 
ac, since by dispersing particles the density of pairs of 
neighboring particles (wf ) decreases whereas the density 
of active sites {ujf) increases. Our results show that for 
sufficiently rapid diffusion, D ^ \, the critical value of 
ac increases without limit, that is, ac — > 00, in agreement 
with the results obtained by Dickman p^ . 

To determine the exponent /3, we have simulated the 
conserved PAM in the supercritical regime for some val- 
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FIG. 4: log-log plot of A = cic — Q versus the density p, in the 
supercritical regime for the conserved PAM for some values 
of diffusion rate. The straight line has a slope 3.61. 



ues of the diffusion rate D. In this regime, the density 
p ^ n/L \s kept fixed for a size system L. We have used 
a lattice with L — 10000 sites and a number of Monte 
Carlo steps ranging from 10^ to 10^. The exponent j3 
is obtained from the log-log plot of IS. = ac ~ a versus 
p as shown in the Fig. |4l We used the values of ac 
calculated by using the equation (fTTjl . For 73 = 0, we 
obtained 1//? = 3.61(3) that is in excellent agreement 
with the value 1//3 = 3.61(1) obtained for the CCP [|]. 
The simulations performed at D = 0.1 and D = 0.3 give 
exponents also compatible with this value. 



The competition between diffusion and reaction pro- 
cess may bring great changes in the phase diagram. For 
example, in the multiple-creation contact processes, the 
transition becomes first order for high enough values of 
diffusion [1, [1]. In the PCP, several works have reported 
a change in the universality class for nonzero values of 
diffusion rate. For the ordinary TAM, previous works 
[To, 11] have shown that in the presence of diffusion, the 
system exhibits a reentrant phase diagram as can be seen 
in Fig. ini 

In references [13, [Hi, the processes of diffusion, cre- 
ation and annihilation are chosen with probabilities Z?*, 
(1 - D*)\/{\ + 1) and (1 - D*)/{\ + 1), respectively. 
To compare ordinary results with ours, it is necessary to 
convert the parameters. The relation between D* and 
A and a are given by 

a 

and 



(16) 



D' 



(17) 



For values of the diffusion rate D* higher than -D*jax — 
0.587, there is no phase transition and the system dis- 
plays only active state. For D* < D^^^, an active state 
is also possible for very small values of A. This happens 
because there are few isolated particles. Since the system 
annihilates only triplets of neighboring particles, these 
sparse particles are able to "survive". On other hand, 
the probability that a new particle is created is very low 
because A is very small. 



B. Constant particle number ensemble 



V. TRIPLET-ANNIHILATION CONTACT 
MODEL 

A. Constant rate ensemble 

The ordinary TAM is composed of the creation of par- 
ticles and annihilation of three neighboring particles. It 
is represented by the chemical reactions: 

+ A^A + A, (13) 

A + A + + + (14) 

describing the catalytic creation and annihilation of a 
triplet of particles, respectively. The quantity cof is also 
given by Eq. ([5]) and the rate ojf is given by 

= Vi-iViVi+i- (15) 

Without diffusion, a continuous phase transition to an 
absorbing state occurs at ac = 0.1488(2) [T^, [HI with 
critical exponent /? = 0.2765. 



The Monte Carlo simulation of the conserved TAM 
is performed as follows: The jumping process and the 
diffusive process are chosen with probabilities 1 — D and 
D, respectively. If the jumping process is chosen, then 
a triplet of neighboring particles is chosen at random. 
Next, we choose three empty sites each one surrounded 
by at least one particle. The three particles belonging 
to the selected triplet jump to these three active sites. 
If the diffusive process is chosen, then a particle chosen 
randomly is moved to one of its neighboring site, provided 
it is empty. The quantities tof and are calculated by 
using the Eqs. ^ and l[T5|) . The rate a is evaluated 
using formula Q with £ = 3. 

To determine the critical value ac we have simulated 
the conserved TAM in the subcritical regime. We assume 
a behavior given by Eq. (jlip for large particle number 
n. For D = 0, the values of a versus n are shown in 
Fig. [5l The linear extrapolation gives the critical rate 
ac — 0.14898(5), which agrees very well with the value 
ac — 0.1488(2) obtained for the ordinary version [lo| . 

For the values of diffusion used here, we have also ob- 
tained values of I /dp consistent with 1.338(6), as shown 
in the Fig. H 
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FIG. 5: Values of a versus the number of particles n for the 
conserved TAM in the absence of diffusion for an infinite sys- 
tem. The line corresponds to the linear extrapolation using 

Eq. unj- 




FIG. 6: Phase diagram for the conserved TAM (full circles) 
in D* versus X^^^ space. The supercritical and subcritical 
regimes are separated by a critical line.. For comparison, we 
show the results obtained by Dickman (squares). 
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FIG. 7: log-log plot of A = — a versus the density p, in the 
supercritical regime for the conserved TAM for some values 
of diffusion rate. The straight line has a slope 3.61. 

priate to determine the transition line. In the ordinary 
ensemble, isolated particles create new particles leading 
to the appearance of triplet even at very small rates and 
small number of particles. In the conserved ensemble, on 
the other hand, it is necessary the existence of at least 
one triplet for the occurrence of a creation-annihilation 
process. In the absence of three adjacent particles, only 
the diffusive process will produce a triplet. If the particles 
are scattered the reunion of three particles may not oc- 
cur. In this case the determination of a will be affected 
by low statistics. We remark that the reentrant phase 
diagram for the TAM was confirmed by time-dependent 
numerical simulations [ll|, which critical exponents be- 
long to the DP universality class. 

In Fig. [7]we show the log-log plot of A = Ofc — a versus 
p for some values of diffusion rate in the supercritical 
regime using a lattice of size L = 10000 and from 10^ 
to lO*" Monte Carlo steps to evaluate the averages. For 
D = 0, we obtained 1//3 = 3.60(3) that is in agreement 
with the value 1//9 = 3.61(1) obtained for the CCP (3]. 
The values for 1 //3 obtained from simulations a,t D — 0.05 
and 13 = 0.1 are also compatible with this value. 



Using the same procedure for nonzero values of D we 
have built the phase diagram, shown in Fig[Bl We used 
the Eqs. ([5]) and (fTT)) to convert the rates. The results 
obtained by Dickman l3,[IH are also plotted in the same 
figure for comparison. 

Our results are in good agreement with those obtained 
for the ordinary ensemble. For example, the results ob- 
tained by Dickman pT| for the maximum critical diffu- 
sion DJ^ax is 0.587 and its respective critical creation rate 
Ac is 0.1. Our results for this point of the phase diagram 
are i'max = 0.589 and Ac = 0.11. Concerning the reen- 
trant phase, the conserved ensemble seems to be inappro- 



VI. PAIR CONTACT PROCESS 

A. Constant Rate ensemble 

The ordinary PCP [T^ is a nonequilibrium model 
which, like the contact process (CP), exhibits a phase 
transition to an absorbing state, but differently from this 
one, the PCP possesses infinitely absorbing states. Nu- 
merical and theoretical studies indicate that the PCP 
(without diffusion) also belongs to the DP universality 
class. 
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The PCP is represented by the chemical reactions: 

+ A + A^A + A + A, (18) 



A + A^ + 0, 



(19) 



P 



a 



a + 1 

B. Constant particle number ensemble 



(20) 



0.0849 



0.0847 - 



0.0845 - 



describing the catalytic creation of a particle and anni- CX 
hilation of a pair of particles, respectively. Notice that it 
is necessary two particles to create a new particle. 

The PCP is defined by the following rules. A pair of 
neighboring particles is chosen at random. With proba- 
bility p it is annihilated and with probability 1 — p a new 
particle is created in one of its nearest neighbor sites. 
The dynamics is governed by pairs of neighboring parti- 
cles, instead of isolated particles. As a consequence, any 
configuration absent of pairs of neighboring particles is 
absorbing. The order parameter is the density of pairs of 
neighboring particles, instead of the density of particles. 

Several studies [l2, [H, H) [13 show that the one di- 
mensional PCP exhibits a second order transition to an 
absorbing state at pc = 0.077090(5) llf\. Close to the 
critical point the density of pairs of neighboring parti- 
cles pp follows the power law behavior with an exponent 
(3 = 0.2765 li*!. The exponents S, rj and z were also 
obtained for the PCP _13J and they are compatible with 
the typical values of the DP universality class. 

Our aim here consists in showing that the PCP can be 
also described by a dynamics that conserves the particle InA 
number. To compare our results with ordinary versions, 
we should note that the strengths of creation and anni- 
hilation rates are kc and ka are proportional to 1 — p and 
p, respectively. Therefore a = ka/kc = p/{\—p) from 
which follows the relation between the parameters p and 
a used here 
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FIG. 8: Values of a versus the number of pairs of neighboring 
particles rip for the conserved PCP in the absence of diffusion 
for an infinite system. The line corresponds to the linear 
extrapolation using Eq. 




The Monte Carlo simulation of the conserved PCP 
without diffusion is performed as follows: A pair of neigh- 
boring particles is selected at random. Next, two active 
sites are chosen at random. Here, active sites are empty 
sites surrounded by at least a pair of neighboring par- 
ticles at the same side. The particles belonging to the 
selected pair jump to the two chosen active sites. The 
rate a is evaluated using Eq. ^ where the creation rate 
is given by 

= (1 - ni)-\\ (21) 

and the annihilation rate uj^ is given by the expression 

To compare the results coming from distinct ensem- 
bles, we have also simulated the ordinary pair contact 
process in one dimension. For both ensembles, we used a 



FIG. 9: Log-log plot of A = Qc — a versus the density of pairs 
of neighboring particles pp, in the supercritical regime for the 
conserved PCP in the absence of diffusion. The straight line 
has a slope 3.61. 



lattice of size L = 10000 and from 10^ to 10^ Monte Carlo 
steps to evaluate the averages. For example, simulating 
the ordinary version for the value of rate a — 0.01 we ob- 
tained the mean density of particles p — 0.96902(1) and 
the mean density of pairs of neighboring particles pp = 
0.94920(2). In the conserved version for p = 0.96900 we 
obtain the averages Pp = 0.94918(2) and a = 0.0100(1). 
Small discrepancies are due the fact that the system size 
is finite. 

In absence of diffusion, any configuration without pairs 
of neighboring particles is absorbing and, therefore the 
conserved PCP also possesses infinitely absorbing states. 
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To determine the critical rate ac we have studied the 
system in the subcritical regime. Since a configuration 
absent of pairs of neighboring particles is absorbing, dur- 
ing the simulations the system may fall into an absorbing 
state. Whenever this happens, we allow isolated parti- 
cles to jump to empty sites surrounded by one particle, in 
order to "create" pairs of neighboring particles. In Fig. 
|5] we show the number of pairs of neighboring particles 
and its respective value of a in the subcritical regime. 
An linear extrapolation in n^^ gives ac = 0.08353(5). 
Using Eq. ^ we find = 0.07709(5) that is in excel- 
lent agreement with the value Pc = 0.077090(5) for the 
ordinary pair contact process [l7| . 

To calculate the fractal dimension dp for the conserved 
PCP, we adopted the following procedure. We measured 
the maximum distance R between two pairs of neighbor- 
ing particles of a subsystem of fixed size L. The log- log 
plot of R versus Up is shown in Fig. [51 The value of 
the inverse of the fractal dimension is 1/dp — 1.33(1), in 
very good agreement with the value 1.338(6) 0. 

Fig. [S] shows a log-log plot of A = ac — a versus the 
density pp. The slope of the straight lines fitted to the 
data points for the conserved PCP has slope 3.61(3), in 
good agreement with the value 3.61(1), obtained for the 
CCP. 



VII. CONCLUSION 

We have analyzed, by numerical simulations, three 
one-dimensional nonequilibrium models in which parti- 
cles belonging to a cluster of £ particles jump to £ dis- 
tinct active sites. They are conserved versions of models 
originally defined in the constant rate ensemble. Our ap- 
proach is general in the sense that any process in which 
the process of creation and annihilation of particles are 
mutually exclusive can be described by a dynamics that 
conserves the number of particles. Our results for the 
three models studied here show that not only universal 
quantities but also nonuniversal parameters are in excel- 
lent agreement with their respective ordinary versions. 
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